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Abstract 

The free energy and the equation of state of isospin- asymmetric nuclear matter are cal- 
culated at finite temperature up to three loop order in the framework of in-medium chiral 
perturbation theory, systematically incorporating one- and two-pion exchange dynamics to 
this order. Effects from the 27r-exchange with explicit A-isobar excitation are included, as 
well as three-body forces. We construct the phase diagram of nuclear matter for differ- 
ent proton fractions Xp and investigate the dependence of nuclear matter properties on the 
isospin- asymmetry. A detailed study of the liquid-gas phase transition is performed. For 
isospin- symmetric nuclear matter we find a critical temperature of 15.1 MeV; as the isospin- 
asymmetry is increased, the liquid-gas coexistence region decreases until it disappears at 
Xp ^ 0.05, while nuclear matter becomes unbound at Xp ^ 0.12. The quadratic expansion of 
the free energy in the asymmetry parameter ^ is a good approximation at low temperature 
even for large S. An estimate of chiral four-body correlations in nuclear and neutron matter 
is also performed and the corrections from such four-body interactions are found to be small. 

Keywords: chiral perturbation theory; nuclear matter; equation of state; thermodynamics; 
phase diagram; liquid-gas phase transition. 

1 Introduction and framework 

Theories and models that describe the thermodynamic properties of nuclear matter play an 
important role in applications to heavy-ion collisions and astrophysics. The plateau exhibited by 
the caloric curve of the nuclear fragments in nucleus-nucleus collisions is interpreted as the trace of 
a first-order liquid-gas phase transition [HIlllS]. Concerning astrophysics, the recent observation 
of a two-solar-mass neutron star [4] provides important constraints for the equation of state of 
strongly interacting matter, as it would rule out many exotic models which tend to produce too 
soft equations of state. 

Within the last two decades a novel approach to the nuclear many-body problem has been 
developed, based on the understanding of the nucleon-nucleon interaction in the framework of 
chiral effective field theories O El [71 [8]. This approach relies on the scale separation between 
the short-range physics and the long- and intermediate-distance dynamics, the latter responsible 
for the attractive part of the nucleon-nucleon interaction described primarily by 27r-exchange 
dynamics. While the long- and intermediate-range, pion-driven dynamics (in the presence of the 
nuclear medium) are treated explicitly, the unresolved short-distance physics is encoded in a few 
contact terms with coefficients fixed in order to reproduce selected known bulk properties of nuclear 
matter. Many other ground state and single-particle properties and nuclear thermodynamics 
emerge as predictions. 

In the following we briefly recall the theoretical framework of in-medium Chiral Perturbation 
Theory (ChPT) developed so far for isospin-symmetric nuclear matter. In the present work we 
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Figure 1: (a): l7r-exchange Fock diagram; (b): iterated l7r-exchange Hartree diagram; (c): iterated 
iTT-exchange Fock diagram. Diagram (a) includes two medium insertions, diagrams (b) and (c) 
involve two and three medium insertions. 

generalize such calculations to isospin-asymmetric nuclear matter with an arbitrary proton-to- 
neutron ratio and perform a systematic analysis of the thermodynamic properties of the resulting 
equation of state. This short summary of the basic framework will be restricted to the essentials 
since the more detailed discussion can be found in Refs. [6^, [7l[3|- 

The relevant active degrees of freedom for the description of the nuclear many-body problem 
are nucleons and pions. In- medium ChPT is applicable to nuclear many-body systems as long 
as the Fermi momentum, kp^ is small compared to the scale of spontaneous chiral symmetry 
breaking in QCD, ~ "^tt/t^ ~ 1 GeV, where /^^ is the pion decay constant in vacuum. The 
characteristic range of nuclear momenta is comparable to about twice the pion mass {kpo ^ 2m-j^^ 
where kpo 263 MeV is the Fermi momentum at the saturation density of nuclear matter, 
po — 0.16 fm~^). This implies that Itt- and 27r-exchange processes have to be treated explicitly in 
the presence of the nuclear medium. Moreover, to complete the picture it is mandatory to include 
the A(1232)-isobar excitation as an explicit degree of freedom, because the splitting between the 
A-isobar mass and the nucleon mass, A = 293 MeV, is again comparable to 2m-j^. Consequently, 
the "small scales" that appear in our calculation scheme are the Fermi momentum of the nucleons, 
kp, the pion mass m^^, and the A (1232) — N mass splitting. The energy per particle at T = is 
then derived as a systematic expansion in powers of kp- The expansion coefficients are non-trivial 
functions of the ratios /ci^/m^^ and of A/m^^. 

Explicit Itt- and 27r-exchange dynamics in the medium are so far treated up to three-loop order 
in the energy density [6, 7, 8 . The basic ingredient in performing calculations at finite baryon 
density is the in-medium nucleon propagator. Consider first the case of zero temperature. Instead 
of the usual decomposition of the propagator into a particle and a hole part, it is convenient to 
use the following equivalent representation: 



where = (po^p) and Mn = 939 MeV is (free) the nucleon mass. The first term is the vacuum 
propagator of the nucleon, while the second term, the medium insertion, takes into account the 
presence of a filled Fermi sea of nucleons. The calculation is then organized according to the 
number of medium insertions in the diagrams representing the energy density. Diagrams with no 
medium insertion give an unobservable shift of the vacuum energy. Diagrams with one medium 
insertion renormalize the nucleon mass to its physical value and provide a description at the level 
of a (non-interacting) Fermi gas. Diagrams with two and three medium insertions produce the 
interesting many-body effects. Two-body terms involving the nucleon- nucleon T- matrix come 
from diagrams with two medium insertions. Three medium insertions generate three-body terms 
such as Pauli-blocking effects on the two-body terms, and so forth. Convergence in the series of 
medium insertions, or powers of /c^, at the three- loop level is realized for /ci? <C ~ 1 GeV as 
long as four- nucleon correlations are not prominent. Effects of four-body correlations are still an 
open issue to be explored in further (ongoing) studies [10 . Genuine chiral four-body contributions 
to the energy per particle are examined later in Section |5j They turn out to be negligibly small. 
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Figure 2: Examples of irreducible 27r-excliange diagrams to be evaluated with two medium inser- 
tions. 
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Figure 3: 27r-excliange diagrams with single and double A(1232)-isobar excitation. 



Fig. [T] shows the Itt- and iterated l7r-exchange diagrams evaluated with up to three medium 
insertions. The iterated one-pion exchange terms include the second-order tensor interaction 
that is well known to be very important in nuclear systems. Fig. |2] displays irreducible 27r- 
exchange diagrams contributing to the energy density. The detailed discussion and terminology 
concerning iterated l7r-exchange and irreducible 27r-exchange diagrams in the context of ChPT, 
both in vacuum and in- medium, can be found in Refs. [6l[7l[8]. The 27r-exchange diagrams with 
inclusion of A-isobar excitations as intermediate states are shown in Fig. |3] Medium insertions 
on the second and third of these diagrams, cutting the intermediate nucleon line, lead to the 
three-nucleon term of Fig. [4j 

Unresolved short-distance physics is encoded in two contact terms, a momentum independent 
term with coefficient Bs and a p^-dependent term with coefficient [8 . For a good description of 
nuclear matter properties an additional attractive three-nucleon contact interaction (controlled by 
the parameter see Eq.(24)) is required, which partly counteracts the strongly repulsive p^-term 
in the energy per nucleon from the three-body diagram with A-isobar excitation. After tuning the 
minimum of the energy per particle of symmetric nuclear matter to —16 MeV, one can fix B^ = 0. 

The next step is the extension to finite temperatures. The ordering scheme in the number of 
medium insertions is now transcribed to the free energy density as a function of density p and 
temperature T. The free energy density is written as the sum of convolution integrals: 

oo oo oo 

pF(p,T)=4 1 dppJClip) d{p)^ J dpi I dp2 IC2{puP2)d{pi)d{p2) 



oo cx) oo 

+ y dpi y dp2 J dp3lC3{pi,p2,P3)d{pi)d{p2)dips) + pAip,T) , (2) 







Figure 4: Three-nucleon term with intermediate A(1232)-excitation. 
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where T) is the free energy per nucleon and 



exp — 



(3) 



is the density of nucleon states in momentum space. Complete expressions of the kernels Kj can 
be found in|Aj The term ^(p, T), the so-called anomalous contribution, is associated with a Fock 
term involving second order l7r-exchange [9l . This contribution vanishes at T = and has its origin 
in the smoothing of the Fermi surface at T 7^ 0. Its effect is small in the range of temperatures 
considered here. The one-body effective "chemical potential" ji in the distribution d{p) is related 
to the density through 



j dppd{p) . (4) 



p = 4 J dppd{p) . 


The kernel JCi is the contribution from the non-interacting nucleon gas, the kernels JC2 and /C3 
describe the effects of the interactions and include the sum of diagrams in Figs. [l][3] with two- and 
three-medium insertions. The pressure is computed using the standard thermodynamic relation 

P„,T) = ,= f^. (5) 

The calculation is then extended to asymmetric nuclear matter and neutron matter. In com- 
parison with the calculation for symmetric nuclear matter, the only changes required are in isospin 
factors and the introduction of two new contact terms associated with the isospin dependence of 
the short-distance interactions, with coefficients denoted by 5^,3 and 5n,5; they are adjusted in 
order to reproduce the value of 34 MeV for the asymmetry energy at the saturation point [8 . 

Note that for pure neutron matter the Pauli principle forbids a three-body contact interaction. 
The behaviour at low density is dictated by the large neutron- neutron scattering length, a^^ 19 
fm [ll|ll2,, and a resummation of the short-distance NN interaction to all orders in the presence 
of the medium is required for a realistic description at low densities {pn ^ 0.02 fm~^)[13 . The full 
implementation of these effects is currently in progress. 

In the present paper we extend the chiral three-loop calculation of isospin-symmetric nuclear 
matter to the isospin-asymmetric case and present a systematic investigation of the phase diagram 
of the corresponding equation of state for different values of the proton-to-neutron ratio. The Fermi 
seas of protons and neutrons are now filled unequally. The propagator ([T]) at T = is split into 
proton and neutron contributions using isospin-projectors: 

eikp - IpI) ^ 0{k, - IpI) + ^ 0{K - IpI) , (6) 

with kp and kn denoting, respectively, the Fermi momenta of protons and neutrons. The "one- 
body" chemical potentials of the two nucleon species in Eq. (|3| are now different. At finite 
temperature we replace 

d{p) — > dp{p) + dn{p) , (7) 

where dp and dn are the proton and neutron distributions. Each diagram now involves the sum 
of all possible combinations of proton and neutron medium insertions with their specific isospin 
factors. In the appendices we list the complete expressions for the diagrammatic expansion of the 
free energy density, both for isospin-symmetric (|A|) and asymmetric matter ([B|). 

This paper is organized as follows. In section |2] we present and discuss the equation of state 
of isospin-symmetric nuclear matter for different temperatures, featuring the liquid-gas phase 
transition. In section [3] the calculations are extended to the isospin-asymmetric case, focusing on 
the behaviour of nuclear matter and its thermodynamic properties with varying proton-neutron 
asymmetry. In section |4] we study the asymmetry free energy, its dependence on density and 
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Figure 5: Free energy per particle of isospin-symmetric nuclear matter as a function of baryon 
density p for different temperatures. The dotted line indicates the non-physical behaviour of the 
free energy in the liquid-gas coexistence region. The physical free energy (solid lines) at low 
temperatures (T < 15 MeV) is obtained using the Maxwell construction. 

temperature and the validity of the parabolic approximation for the free energy as a function of 
the asymmetry parameter, S = {pn — Pp)/{pn + Pp)- Finally, section [t] contains a summary and 
the conclusions. 

2 Equation of state of isospin-symmetric nuclear matter 

We start by recalling and extending the results found for isospin-symmetric nuclear matter 
in Ref. ^ using Eq. ([2|. The free energy per particle T) as a function of density p for 
a sequence of temperatures up to 25 MeV is calculated using Eq. ([2| with the input for the 
interaction kernels ICn specified in[Aj The result is shown in Fig. [5j The dotted lines indicate the 
non-physical behaviour of the equation of state in the liquid-gas first-order transition region. This 
part of each curve is substituted by the physical one (solid lines) using the Maxwell construction. 

At T = the free energy equals the internal energy of the system. The minimum of the curve 
(the saturation point) is located at Eq = —16.0 MeV, po = 0.157 fm~^, kpo = 1.33 fm~^ = 262 
MeV. At finite temperatures the free energy displays a singular behaviour for p ^ 0; this is a 
well-known generic feature that the present calculation shares with other types of many-body 
calculations [15] . 

The energy per particle of symmetric nuclear matter is usually expanded around the saturation 
point: 



with K the compression modulus. From the curve at T = we extract the value K 300 
MeV, slightly larger than the values deduced from relativistic mean field models [25 and from the 
systematics of nuclear monopole resonances [16]. 

Fig. [6] shows isotherms of the pressure P(p,T) (see Eq. ^) as a function of density p. The 
emerging picture is qualitatively reminiscent of a van der Waals gas with its generic liquid-gas 
first-order phase transition. Indeed, chiral A^A^ dynamics generates intermediate range attractive 
interactions, such as 27r-exchange with intermediate A-excitation, that display a characteristic 




(8) 
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Figure 6: Pressure isotherms as a function of density for symmetric nuclear matter displaying a 
first-order liquid-gas phase transition. The dotted lines at low temperature show the non-physical 
behaviour of the isotherms in the transition region. The physical pressure is calculated using the 
Maxwell construction. The dashed line delimits the boundary of the coexistence region. The dot 
indicates the critical pint (Tc :^ 15.1 MeV, pc — Po/3). 



^-2m^r bchaviour at intermediate (1-2 fm) distances (see Ref. [IT]). Such mechanisms account 
already for roughly half of the attraction required to bind nuclear matter at T = 0. The liquid-gas 
phase transition is then the result of a sensitive balance between intermediate range attraction 
and short-range repulsion, the latter represented by contact terms encoding (non-perturbative) 
physics not resolved in detail at the range of Fermi momenta hp relevant to the present study. 

The critical temperature for the liquid-gas phase transition is found at Tc — 15.1 MeV in 
this calculation. For T < Tc the usual Maxwell construction is applied, keeping the pressure 
constant in the liquid-gas coexistence region. Empirical values of this transition temperature 
deduced from multifragmentation and fission measurements locate Tc between 15 and 20 MeV 
[19]. Phenomenological Skyrme model studies [20 gave a similar range for Tc. 

In order to assess possible uncertainties of the thermodynamic properties in our calculation we 
have considered variations of the parameters Br^ and Within the ranges —0.5 < < 0.5 and 
—0.8 < C ^ —0.7 the critical temperature varies from Tc = 14.6 MeV to Tc = 15.8 MeV, i.e. this 
variation implies at most a change of ±5%. 

The discussion of the equation of state is completed by displaying the pressure as a function 
of the nucleon chemical potential (including the free nucleon mass), 

P = Mn^(^I^P^^F{p,T) , (9) 

m Fig.|7l The non-physical (dotted) curves in P(p, T) of Fig. |6] find their correspondence in the 
double- valued behaviour of P as function of p at temperatures below Tc. In this coexistence region 
the actual pressure and chemical potential are constant and given by the intersection point of the 

■"^We briefly remind the main steps of the Maxwell construction. We plot flrst the free energy per particle F(p) 
as a function of the volume per particle v = 1/p. The region with negative curvature is unphysical and is eliminated 
through the double tangent construction. We search for the two points of the curve with the same tangent. This 
tangent line is the physical free energy in the phase coexistence region and its negative slope gives the physical 
pressure. 
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Figure 7: Pressure isotherms as functions of the nucleon chemical potential /i for isospin-symmetric 
nuclear matter. The dotted double- valued region of the curves at temperatures below Tc ^ 
15.1 MeV corresponds to the non-physical behaviour of the equation of state in the liquid-gas 
coexistence region. In this region the actual pressure and chemical potential are constant and 
determined by the Maxwell construction. The dot indicates the critical point. 



double-valued loop according to the Maxwell construction. The temperature at which this loop 
reduces to a point is the critical temperature Tc. 

All relevant informations about the phase transition are collected in the T — /i, P — T and T — p 
phase diagrams (Fig. [8|. In Fig. 8a we display temperature versus chemical potential (T — /i) and 
in Fig. 8b pressure versus temperature (P — T). In these diagrams the coexistence region gets 
projected onto a first-order phase transition line. The first-order transition region terminates at 
the critical point indicated by the dot. The following critical values of thermodynamic quantities 
(pressure, baryon chemical potential and density) are found: ^ 0.261 MeV fm"^ fic - 914.7 



MeV and pc ^ 0.053 fm ^ at T^ = 15.1 MeV. The T-p phase diagram, Fig. [8c) clearly illustrates 
the extension of the gas-liquid coexistence region. 

At T = the third law of thermodynamics gives important constraints about the slope of 
the boundaries of the transition region. It implies that the boundaries of the transition region 
at T = have an infinite slope in the T — p and T — p diagrams and a zero slope in the P — T 
diagram. Moreover, the chemical potential at T = is known and given by the total energy per 
particle at the saturation point, i.e. p = Mn + — 923 MeV. 



3 Equation of state of isospin-asymmetric nuclear matter 

In this section we investigate the dependence of nuclear matter properties on the isospin- 
asymmetry. The amount of asymmetry is given in terms of the proton fraction, Xp = Z/A = 
Pp/{Pp + Pn)- The calculations are again based on the expansion of the free energy density, 
Eq. ([2|, with input specified in[B) 

The resulting equations of state for different values of the proton fraction, Xp = 0.3, 0.1 and 
0, are shown in Fig. [9j The limiting case Xp = corresponds to pure neutron matter. The left 
column of Fig. [9] displays the free energy per particle, the right column shows the pressure, both 
as a function of the baryon density p = Pn -\- Pp- 

As the neutron-proton asymmetry increases, the free energy and the pressure increase at given 
density, indicating the reduced binding in neutron-rich matter. Pure neutron matter is unbound. 
The systematics displayed in Fig. [9] is almost entirely controlled by the isospin dependence of the 
one- and two-pion exchange forces once the relevant contact terms are fixed to yield the empirical 
asymmetry energy at T = 0, about 34 MeV. The coexistence region of the liquid-gas phase 
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Figure 8: Phase diagrams of symmetric nuclear matter. The dot indicates the critical point. 
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0.26 
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0.053 


0.4 
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0.25 
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0.053 


0.3 
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0.22 


920.1 


0.049 


0.2 


10.9 


0.16 


927.2 


0.044 


0.1 


6.4 


0.08 


937.6 


0.037 


0.053 





0.02 


943.6 


0.032 



Table 1: Critical point as a function of the proton fraction Xp — Z/A: listed are the value of 
temperature, pressure, chemical potential and density. 
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Figure 9: Isospin- asymmetric nuclear matter: free energy per particle (left column) and pressure 
isotherms (right column) as a function of nucleon density for proton fractions Xp = 0.3, 0.1 and 
0. The coexistence region, with boundary delimited by the dashed line, gradually diminishes with 
increasing isospin- asymmetry until it disappears at Xp = 0.053. The dot indicates the critical 
point. 
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Figure 10: Dependence of energy per particle and saturation point of nuclear matter at T = on 
the asymmetry. The solid lines represent the energy per particle as a function of nucleon density 
P = pn Pp for different proton fractions Xp = Z/A. The dashed line shows the trajectory of the 
saturation point as Xp varies. For < 0.12 the energy is always positive. 



transition shrinks with decreasing proton fraction until it disappears and only the (interacting) 
Fermi gas phase remains. 

To complete the picture, consider the evolution of the saturation point, defined as the minimum 
of the energy curve at T = 0, in Fig. 10 Starting from its minimum for symmetric nuclear matter, 
£^0 — —16 MeV at po ^ 0.157 fm~^, the binding energy per nucleon is reduced continuously with 
decreasing Xp until it vanishes at a proton fraction Xp c:^ 0.12. Beyond this point neutron-rich 
matter at T = is unbound. 

Phase diagrams of isospin- asymmetric nuclear matter are shown in Fig. 11 for different values 
of the proton fraction, demonstrating how the matter evolves with increasing asymmetry. The 
dashed line shows the evolution of the critical point and its disappearance at Xp ~ 0.053. At this 
particular proton fraction the coexistence region reduces to a point, as can be seen in the T — p 
diagram, meaning that the liquid-gas phase transition does not take place any more. Neutron-rich 
matter with Xp < 0.053 is always in a gaseous phase. The values of the thermodynamic quantities 
at the different critical points are listed in Tab. [T] Note that the critical point disappears at a 
small but finite pressure. 

The following features of the phase diagram in Fig. 11 are worth noting. Above a proton 
fraction close to Xp 0.1, the gas-liquid coexistence region starts at zero density. The T = 
boundary at which this coexistence turns into a Fermi liquid covers a density range between about 
0.5 Po and normal nuclear matter density po, as one moves from Xp = 0.1 to 0.5. A qualitative 
change in the behaviour of the coexistence region takes place at Xp = 0.12, the point at which the 
binding energy at T = vanishes. Following the dashed line in Fig. 10 that represents absolute 
minima of the energy per particle, one observes that for Xp < 0.12 there is still a local minimum in 
E{p^Xp)^ but the absolute minimum is now located at p = 0. Consequently, neutron-rich nuclear 
matter in the vicinity of Xp < 0.12 is a gas at very low density and T = 0, and then enters the 
coexistence region as the density increases. In the range 0.053 ^ Xp ^ 0.12 nuclear matter is not 
self-bound but it can still have a liquid-gas phase transition. 

The present chiral thermodynamics framework for nuclear matter is, of course, oversimplified 
at low densities where nuclear clustering takes place. A detailed study combining the appearance 
of light (deuteron, triton and helium) clusters with relativist ic mean field phenomenology [22^ 
suggests, however, only modest changes of the T — p phase diagram, such as a shift in the position 
of the critical point with changes of less than 10% in Tc and less than 1% in when cluster 
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Figure 11: Phase diagrams of isospin-asymnietric nuclear matter for different proton fractions. 
The dashed hne shows the evolution of the critical point. 
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Figure 12: Asymmetry free energy per particle as a function of the density for different tempera- 
tures. In the left panel 6 = 0.2 {xp = 0.4), in the right panel S — 0.6 [x^ = 0.2). 



formation is incorporated. 



4 Asymmetry free energy 

Introducing the asymmetry parameter 

S = {pn- pp)/p =l-2xp , (10) 

the free energy per particle can be expanded in powers of S around the free energy of isospin- 
symmetric nuclear matter: 

F{pp, p„, T) = F{p, T) + A{p, T) 5^ + 0{S^) . (11) 

This defines the asymmetry free energy per particle, T). The expansion of F involves only 
even powers of S as long as we ignore isospin-symmetry breaking effects. In this limit nuclear 
matter is invariant under the interchange of protons and neutrons. 



In Fig. 12 we show the behaviour of the asymmetry free energy T) for S = 0.2 {xp = 0.4) 



and S = 0.6 {xp = 0.2) as a function of the nucleon density for different temperatures. We note 
that the asymmetry free energy is sensitive to the temperature only at low densities, p < 0.1 fm~^. 

A test of the validity of the parabolic approximation is shown in Fig. [Tsj where we plot the 
free energy difference with respect to isospin-symmetric nuclear matter as a function of for 



different densities p = pn ^ Pp- At T = (left plot of Fig. 13) the linear dependence on 5 is 
seen to be realized very well even up to large 5. At higher temperature (T = 20 MeV, right plot), 
a slight bending is observed especially at low density. In summary, Eq. ( pTj ) is confirmed to be 
a good approximation of the free energy; the term of order 5^ is generally negligible for most 
applications, even for large values of 5. This feature has also been observed in other calculations 
[23[ |24] . Estimates give a value of the quartic term smaller than 1 MeV at the saturation point 

At the saturation density po we have imposed A{pQ^T = 0) 34.0 MeV in our calculation 
in order to fix the contact terms associated with the isospin-dependent part of the interaction. 
For comparison, a relativistic mean-field model [25 constrained by the properties of selected 
nuclei gives A{pq) = 34 ± 2 MeV. A more recent estimate using the same relativistic mean-field 
phenomenology constrained by giant dipole resonances [26] suggests A{pq) in the range (33 - 
37) MeV. A further estimate in the same paper gives the asymmetry energy at lower density, 
p = 0.1 fm~^, in a window between 21 and 23 MeV, slightly lower than our calculated value 
A(0.1fm~^,T = 0) 23.9 MeV. Previous determinations from extensive fits of nuclide masses 
[IZlEg gave A{pq) between 33 and 37 MeV. 
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Figure 13: Free energy per particle of nuclear matter as a function of the asymmetry parameter 
5'^ for different densities. 



Around the saturation point po^ the asymmetry energy at T = can be expanded in powers 
of p — po as follows: 

A{P) = AM + L ^ + ^ f ^) ' + . . . (12) 

3po 2 \ 3po J 

We extract the coefficients L 90.1 MeV and Kas — 153 MeV. The value of L, in particular, is 
compatible with empirical constraints from isospin diffusion which give L = 88 ± 25 MeV [23 . 

For small asymmetries J, the saturation density gets lowered to po{l — 3LS'^/K) and the 
corresponding compression modulus K{S) is often expressed as an expansion in powers of S: 

K{S) = K^ KrS^ + 0{S^) , Kr = Kas-QL , (13) 

where K is the compressibility of symmetric nuclear matter and K^- is called isobaric compress- 
ibility. Our calculated value is = —388 MeV. The empirical determination suffers from large 
uncertainties. A recent result from measurements of the isotopic dependence of giant monopole 
resonances in even-A isotopes gives = —550 ± 100 MeV [29] . 



5 Chiral four-body contributions 

In this section we study four-body contributions to the nuclear matter equation of state. We 
restrict ourselves to the evaluation of a particularly simple class of four-loop in-medium diagrams 
(at T = 0). The leading long-range four-nucleon interaction is constructed by connecting the four 
nucleon lines through exchanged pions which couple together according to the chiral pion-pion 
interaction [30]. However, since the off-shell chiral 47r- vertex is involved in the process this con- 
tribution alone is representation dependent and thus not unique. It has to be supplemented by 
the four-nucleon interaction generated additionally by the chiral A^A/'37r-vertex, where the three 
pions emitted from one nucleon are absorbed on each of the other three nucleons. The combina- 
tion of contributions from the chiral TTTr-interaction and the chiral A/" A/'37r- vertex is representation 
independent and thus gives rise to unique and physically meaningful results. The (isovector) 
Weinberg- Tomozawa NNtttt involves the energies of the exchanged pions which in the present 
application are differences of nucleon kinetic energies. The corresponding four-nucleon interac- 
tion is then a relativistic 1/M^-correction. Note that we do not follow the method of unitary 
transformations of Ref. [30]. In this scheme three-body and two-body forces and even discon- 
nected diagrams give rise to induced four-body forces. We consider here only diagrams related to 
"genuine" four-nucleon forces. 



Fig. 14 shows the in-medium diagrams one obtains by closing the four nucleon lines to two 



rings. The combinatoric factors of these diagrams are 1/8 and 1/2, respectively. Diagrams with 
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Figure 14: In-medium Hartree diagrams related to the leading-order chiral four-nucleon interac- 
tion. Their combinatoric factors are 1/8 and 1/2. The four- medium insertions are indicated by 
the double-line breaks on the nucleon propagators. 



more rings are trivially zero, due to a vanishing spin- trace. The integrals over the Fermi spheres 
factorize and thus can be readily solved. After adding both diagrams one gets the following 
Hartree contribution from the chiral four-nucleon interaction to the energy per particle of isospin- 
symmetric nuclear matter: 



E{p) 



(AN) 



9g\mlu 



2u arctan 2u 



1 

8^ 



ln(l + 4^x2) 



(14) 



with u = kF/rriT^. Note that after adding both diagrams in Fig. 14 the only remainder of the chiral 
TTTT-interaction is a constant factor — 3m^//^. This feature has ultimately lead to the expression 
with a complete square in Eq. (14). The upper curve in Fig. 15 shows the repulsive four-body 
contribution E{p)^^^^ as a function of the nucleon density p = 2/c|./37r^. One observes that up 
to twice normal nuclear matter density, 2po = 0.32 fm~^, the effects from the genuine long-range 
(pion-induced) four-nucleon correlations stay below 0.1 MeV, and thus are negligibly small. To be 
specific, at saturation density one has E{po)'^^^^ 18keV. 

The chiral four-body contribution in pure neutron matter is also of interest. In this case only 
neutral pions are present which leads to a reduced isospin weight factor. Doing the calculation, 
one finds for the energy per particle En{pn)^^^^ of pure neutron matter the same expression 
as in Eq.(14) with an additional prefactor —1/6 and the substitution u = /Cn/m^^, with kn the 



neutron Fermi momentum. The lower curve in Fig. p~5] shows the attractive four-body contribution 
En{pn)^^^^ as a function of the neutron density pn = /c^/37r^. Note that the increase of the Fermi 
momentum (by a factor 2^/^) has effectively compensated the (isospin) reduction factor 1/6. Again, 
with less than 0.1 MeV in magnitude this long-range four-neutron correlations are negligibly small. 



In addition to the two Hartree diagrams shown in Fig. 14 there are (exchange) Fock diagrams with 
one single closed nucleon ring. We do not enter here into the intricate numerical evaluation of 
these (non-factorizable) four- loop diagrams. But from general considerations one knows that Fock 
terms come with opposite sign and reduced spin- and isospin weight factors, and therefore one 
expects some even smaller and partly compensating effects. 

A more comprehensive study of pion-induced four-nucleon forces including also virtual A(1232)- 
isobar excitations in the context of nuclear matter calculations is an ongoing project. In comparison 



to eq, (14) one expects more significant effects from these mechanisms. First, there exist three-ring 
diagrams with higher spin-isospin weight factors and, secondly, the corresponding contributions to 
the energy per particle do not vanish in the chiral limit 0. Note that Ref. [31] has already 

studied such pion-induced A(1232)-excitation mechanisms for the four-nucleon system. 



6 Comparison to other approaches 

In the recent literature several calculations of the thermodynamic properties of nuclear matter 
using different approaches can be found [32l [33l [34l [35l |36] . In the following we comment on some 
of these, based either on the self-consistent Green's function (SCGF) method [34] or on mean field 
theory [36] . 

As pointed out in Ref. [32], many traditional studies of nuclear matter at finite temperature 
in mean field approximation treat the temperature dependence in a naive way. The temperature 
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density p or [fm ] 

Figure 15: Contribution to the energy per particle arising from the long-range chiral four-nucleon 
interaction. Upper curve: isospin-symmetric nuclear matter, lower curve: pure neutron matter. 

dependence arises mainly from the replacement of the step function momentum distribution at 
zero temperature by the corresponding Fermi-Dirac distribution. In this way the temperature de- 
pendence of the phenomenological interactions accounting for the correlations between nucleons is 
ignored. On the other hand, in many-body calculations based on microscopic approaches, such as 
ours, the nuclear medium generates Pauli blocking effects that are weakened with increasing tem- 
perature. Consequently nuclear matter properties and correlations in the medium are temperature 
dependent in a non-trivial way. 

In the SCGF approach the ladder approximation provides the minimal scheme for the de- 
scription of the effective NN interaction in the medium. The three-body Green's function is 
decomposed into a one-body and a two-body propagator and arranged in an integral equation 
for the in- medium T- matrix that sums the ladder diagrams to all orders. The set of equations 
including the T-matrix, the nucleon self-energy and the single-particle Green's function can then 
be solved self-consistently. In Ref. f3^ the equation of state of isospin-symmetric nuclear matter 
has been calculated using such an approach in combination with the CD-Bonn and Nijmegen NN 
potentials. Contributions of three body forces are included via an effective (density-dependent) 
two-body interaction. The pressure P(p, T, S) is found to be strongly dependent on three body cor- 
relations and the liquid-gas coexistence region gets reduced in size when they are included. With 
the CD-Bonn NN potential the critical temperature results in Tc = 12.5 MeV. When using the Ni- 
jmegen potential Tc reduces to 11.5 MeV. The critical density lies in the range pc — 0.09—0.11 fm~^ 
and the critical pressure is about 0.15 MeV fm~^. The results of these calculations, including the 
very low value of Tc, differs significantly from ours (compare with Tab. [T]). 

In Ref. [36] the properties of nuclear matter have been calculated in the framework of rela- 
tivistic mean field theory with density-dependent meson-nucleon couplings accounting for medium 
modifications. This method has been successful in describing many properties of both nuclear 
matter and finite nuclei. For symmetric nuclear matter a critical temperature of Tc = 13.2 MeV 
is found, while the liquid-gas coexistence region vanishes at Xp 0.07. The emerging picture 
is closer to the one we have obtained in the present paper. The main difference is a somewhat 
smaller phase transition region. We recall that the empirical range for the critical temperature is 
Tc - (15-20) MeV [l9]. 
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7 Summary 



The present work extends previous calculations of nuclear chiral thermodynamics from symmet- 
ric nuclear matter to the isospin-asymmetric case. A virial expansion of the free-energy density 
is treated to three-loop order using in-medium chiral perturbation theory. One- and two-pion 
exchange processes in the nuclear medium are computed explicitly with inclusion of A-isobar ex- 
citation. Short-distance physics is encoded in a few contact terms, the only parameters of the 
approach, fixed to reproduce selected properties of normal nuclear matter at zero temperature. 
The predicted nuclear thermodynamics is then entirely determined by in-medium pion-nucleon 
dynamics with no additional free parameters. Correlations in the nuclear medium include two- 
and three-body forces. 

Within this framework the nuclear equation of state for different temperatures is calculated 
for a series of proton fractions Xp = Z/A in the full range < Xp < 0.5 between pure neutron 
matter and symmetric nuclear matter. The evolution of the first-order liquid-gas phase transition 
is systematically investigated. The critical temperature Tc — 15.1 MeV at Xp — 0.5 falls as Xp 
decreases, and reaches Tc = at 0.053, the point at which the liquid-gas transition disappears. 
The free energy is well approximated at low temperature by the usual quadratic expansion in the 
asymmetry parameter 5 = 1 — 2xp^ even for large 5. 

Part of the purpose of the present work is to set realistic constraints for ongoing discussions of 
the QCD phase diagram in the nuclear terrain of moderate baryon densities (up to about twice 
the density of normal nuclear matter) and temperatures up to about 50 MeV. For example, nu- 
clear chiral thermodynamics with its explicit treatment of pionic degrees of freedom in the nuclear 
medium is the appropriate framework that permits to calculate the density and temperature de- 
pendence of the chiral (quark) condensate beyond the Fermi gas limit [33 [38l [39] , while respecting 
at the same time important nuclear physics constraints. 
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A Contributions to the free energy density: 
isospin-symmetric nuclear matter 

In this appendix we summarize the diagrammatic contributions to the free energy density by 
reporting the explicit expressions of the kernels for the convolution integrals Eq. For further 
details and explanations see Refs. [3 [8]. 

One-body kernel 

2 4 

where jl is the effective one-body chemical potential and = 939 MeV is the free nucleon mass. 

Two-body kernels 
• Contact terms: 

4''\p^,P2) = 247r^53^ + 207r^B5^(p? +^1) , (16) 
with B3 = -7.99 and B5 = 0. 
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• iTT-exchange Fock- diagram: 
35i 



1 



m2 + (pi -p2)^ 

2 / 2 , 2m + (Pl +P2)^ 



m2 + (pi -p2)^ 
2m2pip2(p? -^2)^ 



}■ 



(17) 



• Iterated lyr-exchange Hartree-diagram: 

JCi'^^\pi,P2) = ^ ^^J^Z^'' l{Pi +P2)arctan 



Pi +P2 



+ (P2 — Pi j arctan -m^ m 



8 ml^{pi-p2y 



(18) 



• Iterated l7r-exchange Fock- diagram: 

TT 



327r/4 



(Pl+P2)/2m7r 



2 Pi p2 + 



dx 



1 + 2^2 



[(1 + 8x^ + 8x^) arctan X 



bi-P2|/2m7, 



-(1 + 4x^) arctan 2x]} . (19) 



• Irreducible 27r-exchange Fock-diagram: 



(irrF) 



(Pl,P2) 



2m7, 



M^ + (Pl+P2)^ 4piP2 4piP2 / 2 , 2^ 



X < /iln 



+ (pi -p^y 



(20) 



where ImFc, Iin^c, ImVr, IniVFr are the spectral functions of the isoscalar and isovector 
central and tensor A^A^-amplitudes (see section 3 in Ref. [lij- 



Three-body kernels 

• Iterated lyr-exchange Hartree-diagram: 

P1+P2 

IC3 '{Pl,P2,P3) 



bi-P2| 



\pl-pl^q^^2psq\ 
{ml + ^2)2 - \p2 _p2^q2 _ 2p^q\ 



In 



• Iterated lyr-exchange Fock-diagram: 



^t^\pi^P2,P^) 



n i8p 



4piP3 + {pI-pI - ml) In 



ml + {pi +P3) 
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16/, 

4P2P3 + (P3 -P2 - ^1) In 



m2 + (pi -p^Y 
ml + (p2 +P3)^ ^ 

+ (P2 -P3)^ 



P2+P3 



/ 



q 



b2-P3l 



In + ^1 + !!i In 

|p2-/i| i? |pii?+(pi + m2 -p2)/j| 



(21) 



(22) 
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h=Yq{pl-Pl-^')^ R=\/{nil+pj-plf+iml{pl-h^) . (23) 



• 27r-exchange three body Hartree-diagram with single A-isobar excitation: 

/C3 iPuP.,Ps) = [2p,P. (1 + C) + ^ml + ip,+p.mml + (p, - p.)^] 

-n^^ln4±fl±i^;| , (24) 

including the contact term with = —3/4. 

• 27r-exchange three body Fock-diagram with single A-isobar excitation: 

4 

/Cf ^Hpi,P2,P3) = -T^r^[2X(pi,P3)X(p2,P3) + nPl,P3)nP2,P3)] , (25) 

X(p.,P3) = 2,,p3-^ln4±fl±^, (26) 
r(p„,3) = ^ i^pl - Sml - 3pl) + ^(Pl-Pl + fr + ^<Pi 1, rnlHp^+psr 

4P3 16^3 TO^ + (pi-p3)^ 

(27) 

Anomalous contribution 



pA(p,T) 



m: 



CX) oo 

dpi J dp2 J dp3d{pi)d{p2)[27r'^d{p2) -P2] d{p3 
000 



4p2 rni^[pi-p2)^ 



4p2 



(P3 -P2)^ 



with 



dpi / dp2d{pi) 



d{p2 

P2 



{pi +P2)' 



ml + {pi +^2)^ 

00 

d(p) 



dp - 



P 



{Pl-P2f 



(28) 

(29) 
(30) 



B Contributions to the free energy density: 
isospin-asymmetric nuclear matter 

In isospin-asymmetric nuclear matter the Fermi momenta and the "one-body" chemical po- 
tentials of protons and neutrons are different. As a consequence, Eq. Q has to be modified, by 
taking into account that 

• for each closed diagram we sum over all possible combinations of protons and neutrons, each 
combination multiplied by its own isospin factor; 

• the densities of states dp{q) and dn{q) are different due to different jip and jln- Note that 
the momenta to be integrated over are now denoted by the symbol q instead of p. 
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One-body kernel 

The integration over the proton and neutron distributions have to be performed separately, 
resulting in the following replacement in Eqs. ([2| and (15): 



■I' 



4 / dqqlCi{q)d{q) ^2 I dqq \lC^^\q)dp{q) + ic[''\q)d^{q) 



(31) 



where JC^'^\q) = flp 



3Mn SMfr ' 



Two-body kernels 

• iTT-exchange Fock-diagram. 

Using JC^^-' of Eq. (17) the product d{qi)d{q2) for isospin-symmetric nuclear matter is re- 
placed by: 



d{qi)d{q2) — > - [dp{qi)dp{q2) + c^n(<7i)<^n(<72) + 4 ^^(gijdnfe)] • 



(32) 



• Iterated lyr-exchange Hartree-diagram. 

1 



Using JC2^^^ from Eq. (IT8|), replace in the corresponding integral: 



d{qi)d{q2) 



12 



[dp{qi) dp{q2) ^ dn{qi) dn{q2) + dp{qi) dn{q2)\ 



(33) 



• Iterated l7r-exchange Fock-diagram. 

Using IC^^^^ from Eq. (19), replace in the corresponding integral: 



d{qi)d{q2) — > -[S dp{qi) dn{q2) - dp{qi) dp{q2) - dn{q2) dn{q2)] 



(34) 



• Irreducible 27r-exchange. 



For isoscalar amplitude Vc,t in Eq. (20), replace: 



d{qi) d{q2) — > ^ [dp{qi) dp{q2) + dn{qi) dn{q2)] 



(35) 



For isovector amplitude Wc,t in Eq. (20), replace: 



d{qi)d{q2) — > ^[dp{qi) dp{q2) ^ dn{qi) dn{q2) ^ ^dp{qi)dn{q2)] 



(36) 



• Contact terms. 



Expressions for symmetric nuclear matter using /C2^*^ of Eq. (|16|) are replaced as follows: 



Bs d{qi) d{q2) — > {Bs - Bn^s) dp{qi) dn{q2) + ^^^'^ dp{qi) dp{q2) + dn{qi) dn{q2) , (37) 
^5 d{qi) d{q2) — > {B^ - Bn^b) dp{qi) dn{q2) + ^^n,5 dp{qi) dp{q2) + dn{qi) dn{q2) , (38) 



with Bn 3 = -0.95, Bn 5 = -3.58. 
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Three-body kernels 



• Iterated l7r-exchange Hartree-diagram. 

Using /Cg*^^ from Eq. (21), replace in the corresponding integral: 



d{qi) d{q2) d{q3) 



12 



dp{qi) dp{q2) dp{qs) + dn{qi) dn{q2) dniqs) + dp{qi) dp{q2) dn{q^) 



+ dn{qi) dn{q2) dpiqs) + 4 dn{qi) dp{q2) dniqs) + 4dp((7i) dn{q2) dp{q^) . (39) 



• Iterated lyr-exchange Fock- diagram. 

-3 



Using /C3 from Eq. p2|) , replace in the corresponding integral: 



d{qi) d{q2) d{qs 



2 dp{qi) dn{q2) dp{q^) + 2 dp{qi) dn{q2) dniqs) - dp{qi) dp{q2) dp{qs) 



+ 2 {qi ) dp{q2) dp{qs) ^ 2 dn {qi ) dp {q2 ) dn {qs ) - dn {qi ) dn {q2 ) dn {qs ) • (40) 

• Irreducible 27r-exchange Hartree-diagram with single A-isobar excitation. 
Using JC^^^^ from Eq. (24), replace in the corresponding integral: 



d{qi) d{q2) d{qs) 



1 

12 



dp{qs)^dn{qs) dp{qi) dp{q2)^dn{qi) dn{q2)^^dp{qi) dn{q2) 



(41) 



The three-body contact term proportional to is excluded from this formula, because it 
transforms differently: 



C d{qi) d{q2) d{qs) — > C ^dp{qi) dn{q2) [dp (qs) ^ dniqs)] 



(42) 



• Irreducible 27r-exchange Fock-diagram with single A-isobar excitation. 

Expressions for symmetric nuclear matter using KL^^^ of Eq. (25) are replaced as follows: 



2X{qi,qs)X{q2,qs)d{qi)d{q2)d{qs) — > -X{qi,qs)X{q2,qs) x 



dp{qi) dp{q2) dp{qs) + dp{qi) dn{q2) dp{qs) + dn{qi) dp{q2) dp{qs) + 3 dn{qi) dn{q2) dp{qs) 
^dn{qi) dn{q2) dn{qs)^dn{qi) dp{q2) dn{qs)^dp{qi) dn{q2) dn{qs)^^ dp{qi) dp{q2) dn{qs) , 



y{qi^qs)y{q2^qs)d{qi)d{q2)d{qs) — > -Y {qi , qs)Y {q2 , qs) x 

o 

2 dp{qi) dp{q2) dp{qs) - dp{qi) dn{q2) dp{qs) - dn{qi) dp{q2) dp{qs) + 3dn(<7i) dn{q2) dp{qs 
+2 dn{qi) dn{q2) dn{qs)-dn{qi) dp{q2) dn{qs)-dp{qi) dn{q2) dn{qs)^^ dp{qi) dp{q2) dn{qs) 



(43) 



where X{qi^qs) and Y{qi,qs) are given in Eq. (26) and (27) 



Anomalous term 

The anomalous term is given by the sum of two terms transforming in a different way. 
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• Replace in the triple integral in Eq. (28): 
d{qi) d{q2) [27i'^d{q2) - q2\ d{q^) — > 

^|(ip((7i)dp((73){(ip((72)[27r^<ip(g2) - ^2] + 4 (in(<72) [27r^(in(<72) -^2]} 
+ dn(<7i)c^n(<73){c^n(<72)[27r^(in(<72) - ^2] + 4 ^^((72) [27r^<ip(<72) - ^2]} 

-V4.dp{qi)dn{q^){dp{q2)[27i'^dp{q2) - ^2] + c^n(<72)[27r^(in(g2) -92]}}. (44) 



• The subtraction term in Eq. (28) transforms as: 



2n'i{p,T) 2%^{Pp,T) 2nUpn,T) 

with p = pp -\- pn and 



(45) 



00 

^opAPp,n, T) = -2Mm I dq . (46) 



Using Vl'i^{p^T) defined in Eq. (29), ^'i^^p is obtained replacing in the integral: 

d{qi) d{q2) — > ^ [dp{qi) dp{q2) + 2 dn{qi) dp{q2) . (47) 
Using Vl'i^{p^T) defined in Eq. (|29|, ^i^^^ is obtained replacing in the integral: 
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